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We investigate the electronic structure and properties of atoms exposed to a magnetic quadrupole 
field. The spin-spatial as well as generalized time reversal symmetries are established and shown to 
lead to a two-fold degeneracy of the electronic states in the presence of the field. Low- lying as well 
as highly excited Rydberg states are computed and analyzed for a broad regime of field gradients. 
The delicate interplay between the Coulomb and various magnetic interactions leads to complex 
patterns of the spatial spin polarization of individual excited states. Electromagnetic transitions 
in the quadrupole field are studied in detail thereby providing the selection rules and in particular 
the transition wavelengths and corresponding dipole strengths. The peculiar property that the 
quadrupole magnetic field induces permanent electric dipole moments of the atoms is derived and 
discussed. 
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I. INTRODUCTION 



External fields represent an excellent tool to influence and control the motion of quantum systems and in particular 
to generate new, sometimes very surprising, properties of e.g. atoms and molecules. In the case of static homogeneous 
magnetic fields the interest in the electronic structure of the hydrogen atom reached its maximum in the late 1980's 

. This was motivated on the one hand by the experimental accessibility of highly excited Rydberg states via laser 
spectroscopy as well as the astronomical observation of spectra from magnetic white dwarfs possessing a hydrogen-rich 
atmosphere 2] . On the other hand the development of new numerical techniques allowed to address the regime where 
the competition of the Coulomb and diamagnetic interaction leads to unusual and complex properties and phenomena. 
Besides the vivid interest in magnetized structures the hydrogen atom served as a paradigm of a nonseparable and 
nonintegrable system possessing major impact with respect to the development of several fields such as quantum 
chaos, semiclassics of nonintegrable systems and nonlinear dynamics in general (see 0,0,0,0 and references therein). 

In contrast to the case of a homogeneous magnetic field there exist no investigations on the electronic structure 
and properties of atoms in inhomogeneous or trapping magnetic field configurations. Apart from the fact that this 
problem is of fundamental interest, experimental techniques allow nowadays to manufacture magnetic microtraps on 
atom chips 6]. The latter consists of surface-mounted structures carrying microscopic charges and/or currents which 
create a miniature landscape of inhomogeneous electric and/or magnetic fields. In particular it is possible to build 
tight magnetic traps which involve magnetic field configurations with large field gradients up to B ~ 10 8 ^-. Since 
these traps strongly confine the atoms it is to be expected that the electronic structure of excited or highly excited 
Rydberg states is influenced significantly by the spatially varying field. Quadrupole magnetic fields represent a key 
configuration which is used in many setups of traps. We therefore investigate here the electronic structure of atoms 
exposed to a quadrupole magnetic field (see also ref.Q). 

In detail we proceed as follows. In section [n] we present the electronic Hamiltonian of the atom in the quadrupole 
field including the interaction of the spatial as well as spin degrees of freedom with the external field. Section 
IIIII contains a discussion of the remarkable spin-spatial symmetries of the Hamiltonian including unitary as well as 
antiunitary symmetries and the related constants of motion. As a result we encounter a two-fold degeneracy of 
each eigenstate. Section |^ discusses the eigenstates obeying different symmetries and the resulting consequences for 
expectation values of observables. Section Ivl briefly outlines our computational approach. In section IVTI we present 
the results of our numerical investigations establishing new spectral and other properties for low-lying and highly 



* ilesanov@physi.uni-heidelberg.de 
joerg.schmiedmayer@physi.uni-heidelberg.de 

* Peter.Schmelcher@pci.uni-heidelberg.de 
§ Corresponding author 



2 



excited states for weak and strong gradients. The spin expectation values and in particular the spatial distribution 
of the spin polarization of the excited states are studied in detail. Electromagnetic transitions in the quadrupole field 
including dipole strengths and wavelengths are given too. Whereever appropriate a comparison with the case of a 
homogeneous magnetic field is performed. Finally the peculiar property of magnetic field-induced permanent electric 
dipole moments is derived and investigated in depth. Section IVIII contains the conclusions and outlook as well as 
some discussion of relevant aspects going beyond the present investigation. 



II. THE HAMILTONIAN 



In ultracold atomic physics moving atoms in inhomogeneous magnetic fields are considered to be neutral point-like 
particles which couple only through their total angular momentum to the external field 0, 0- ^ n slowly varying 
magnetic fields it is appropriate and convenient to describe the dynamics of the atoms by means of an adiabatic 
approximation, where one assumes that the magnetic moment is oriented either parallel or antiparallel to the magnetic 
field. In this way the interaction potential becomes proportional to the modulus of the magnetic field. As long as 
the spatial extent of the atom is much less than the typical variations of the field and if one is only interested in the 
center of mass motion of the atom the adiabatic approximation can be applied reliably. 

In the present investigation we are interested in the electronic structure of excited atoms exposed to a quadrupole 
magnetic trap or field. Therefore we have to go beyond the above approximation in the sense that the detailed coupling 
of both the electronic spatial as well as the spin degrees of freedom to the external field is taken into account. The 
majority of todays experiments on ultracold atoms deal with alkali-atoms 6]. Alkali-atoms become rapidly hydrogen- 
like if their single valence electron is excited. The detailed structure of the core becomes increasingly irrelevant with 
increasing degree of excitation. We therefore assume that the outer electron is exposed to the field of a single positive 
point charge. Furthermore we neglect the spin-orbit (LS) coupling and the coupling of the nuclear and electronic 
spin, i.e. the hyperfine interaction. These interactions show an r _3 -dependence, i.e. a rapid decay with increasing r. 
Thus for sufficient highly excited states their contributions are negligible compared to those of the Coulomb-potential 
and the interaction with the external field. Nevertheless one can discuss how the inclusion of spin-orbit coupling 
might modify the properties of the system. For hydrogen in a homogeneous field this relativistic effect changes the 
symmetries of the system crucially: S z and L z are no longer conserved separately and only J z remains a conserved 
quantity. The result is a restructuring of the energy spectrum since previously separated symmetry sub-spaces become 
coupled. In particular this manifests itself in the emergence of avoided crossing Q. For the atom in the quadrupole 
field such dramatic effects are not expected since the additional LS'-coupling term does not reduce the symmetry of 
the system. Additionally, considering its r _3 -dependence, the spin-orbit interaction can certainly be treated by means 
of perturbation theory in the regime of spatially extensive and energetically highly excited Rydberg states. 

In the presence of a magnetic field the motion of the center of mass (CM) and internal (electronic) degrees of freedom 
of an atom do not decouple. This holds in particular for the case of a homogeneous magnetic field |?J [HI EI ^3 
and therefore has also to be expected for an inhomogeneous field. However, CM-motional effects on the electronic 
structure become only significant in certain parameter and/or energetic regimes I^ElEl- I n order to approximately 
decouple the CM and electronic dynamics we take advantage of the heavy atomic mass compared to the electron 
mass (ttia ^> m e ). Additionally we assume an ultracold atom whose CM motion takes place on much larger time 
scales than the electron dynamics, even for highly excited electronic states. The Hamiltonian describing the electronic 
motion then reads 
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Where we have assumed that the atomic core (nucleus) is fixed at the trap center which is the origin of the coordinate 
system. The magnetic quadrupole field is given as 

B(f)=b^ y j. (2) 

The gradient b is the only parameter characterizing the steepness of the quadrupole trap. The vector field is rotationally 
symmetric around the z-axis and invariant under the z-parity operation. Figure Q shows intersections of this vector 
field for z = and x = 0, correspondingly. A vector potential belonging to the quadrupole field reads 

A(f) = i B(r) x r = b \ -xz ) , (3) 
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FIG. 1: Vectorial plots of the magnetic quadrupole field J2J). The intersections are the (x, y)-plane for 2; = (a) and the 
(y, z)-plane for x = (b). 



which is a Coulomb gauge (VA(r) =0). In Cartesian coordinates and adopting atomic units 30] (g s = 2) the spinor 
Hamiltonian becomes 



, 5 1 5 ? ; - 6zL z + ^z 2 (x 2 + y 2 ) + I (<t x z + <j y y - 2<r*z) 



(4) 



with L z = —i (yd x — xd y ) and cr x , a y , a z being the Pauli spin matrices (S = |<r). The first two terms of He in equation 
J3J) represent the Hamiltonian of the hydrogen atom without external fields. The third term is the paramagnetic (oc b) 
or Zeeman-term. In contrast to the situation of the atom in a homogeneous field this term now depends not only on 
L z but also linearly on the z-coordinate. The fourth (diamagnetic) term (ex b 2 ) represents a quartic oscillator coupling 
term between the cylindrical coordinates p = y ' x 2 + y 2 and z. In a homogeneous field the diamagnetic interaction 
is a pure harmonic oscillator term proportional to p 2 and yields a confinement perpendicular to the magnetic field. 
Finally the reader should note that the coupling of the electronic spin to the quadrupole field depends linearly on the 
Cartesian coordinates. The latter prevents the factorization of the motions in the spin and spatial degrees of freedom 
and renders the corresponding Schrodinger equation into a spinor equation. This is again in contrast to the situation 
of a homogeneous magnetic field where the spin component along the field is a conserved quantity. 



III. SYMMETRIES 



In this section we investigate the unitary as well as the anti-unitary symmetries and conserved quantities of the 
Hamiltonian (|3J). We will then construct the corresponding eigenfunctions respecting these symmetries. This will 
help us (in Section to derive properties of expectation values. Some basic definitions and features of probability 
densities are also introduced being utilized later to analyze excited atomic states in the quadrupole field. 

For the following discussion it is useful to transform the Hamiltonian (@J) into spherical coordinates. It reads 

H s = -^A r> 0,0- ^ + y^ 4 cos2 sin2 6 +\ r sin OK-br cos 6 (L z + a z ) (5) 
with K being the matrix 

*=(£ e 7)- ^ 

As a consequence of the rotational invariance of the quadrupole field the z-component of the total angular momentum 
J z = L z + S z is conserved, i.e. we have [Hs,J z ] = 0. Additionally we have the discrete symmetry represented by 
the unitary operator P^OP^, i.e. \H s ^P$OP z \ = 0. Here O = a x exchanges the components of a ^-spinor, P z is the 
z-parity operator P z : z — > — z and P^ represents the 'reflection' P^ : <\> — > 2ir — <p. Apart from these symmetries the 
Hamiltonian possesses two generalized anti-unitary time reversal symmetries namely TOP z and TP^. Both of them 
involve the conventional time reversal operator T (T 2 = 1), which corresponds to complex conjugation. Applying the 
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two anti-unitary operators TOP z and TP^ results in the unitary operation P^OP z . The operators TOP z , TP^ and 
P^OPz form an invariant Abelian sub-group. Together with J z they obey the following (anti-)commutation rules: 

[J z , TP^} = {J z , TOP z } = {J z ,P^OP z } = (7) 
[TP h TOP z ] = [TP , P^OP z \ = [P^OP z , TOP z ] = (8) 

Apparently the spin-spatial symmetry operations form a non- Abelian symmetry group. This group generated by 
P(f)OP z and J z is isomorphic to C s . 

In the following we will show, that the interplay of the above symmetries leads to a degeneracy of the eigenvalues of 
the Hamiltonian Q. We consider a state \E, m) which is an eigenstate of the Hamiltonian JJJ with the energy E and 
to J z with the half-integer quantum number m (see below). Since TOP z commutes with Hs the state TOP z \E,m) 
is also an energy eigenstate with the energy E. By letting J z act on this state one obtains 

J z TOP z \E,m) = -TOP z J z \E,m) = -m TOP z \E,m) . 

Thus the state TOP z \E,m) can be identified with \E, —m). Hence the states with the eigenvalues m and — m are 
degenerate. This two- fold degeneracy of each energy level in the presence of the inhomogeneous magnetic field is a 
remarkable feature. It is reminiscent of the Kramers degeneracy of spin - systems in the absence of external fields 

M- 

In the case of an atom in a homogeneous magnetic field the operators L z , P z as well as parity P and Ta z P$ form 
the corresponding set of spatial and time reversal symmetries if the field is oriented along the z-axis. Since these 
operators commute they form an Abelian symmetry group implying the well-known fact that there are no energy level 
degeneracies in the homogeneous field. 

In order to finalize this section we want to remark that most of the above considerations hold not only for the case 
of a fixed nucleus in the trap center but also for the full two-body problem. If the dynamics of both the nucleus and 
the electron is considered the conservation of TP^, TOP z and J z of the two-body system holds. 



IV. J z -, TOP z - AND TP^-EIGENSTATES 

Since the J z operator provides the eigenvalue equation J z \m) = m \m) with half-integer m the spatial representation 
of the spinor eigenfunctions \m) are 



i(ra-f )0 



m) = [ : i(m+i)<b ) ■ (9) 



The eigenfunction \m) can be used to reduce the electron dynamics to a given m-subspace. The quantized motion of 
the electron is then described by the Hamiltonian H m = (m\H \m). We will exploit this fact later when solving the 
stationary Schrodinger equation. 

Eigenstates to the TOP^-operator can be constructed by superimposing two degenerate J^-eigenstates. They read 

\E, ±) TOP * = -J= [\E, m) ± TOP z \E, m)} . (10) 

The corresponding eigenvalue relation is 

TOP z \E,±)=±\E,±). (11) 
Eigenfunctions of the anti-unitary operator TP^ can be constructed analogously: 

\E,±) TP * = -^=[\E,m)±TP^ \E,m}} (12) 

The J z - and TOP^-eigenstates exhibit certain symmetry properties. Exploiting such symmetries often turns out 
to be extremely powerful and can save a lot of computational effort. It is especially useful when one is interested in 
computing matrix elements such as expectation values of certain observables. 

The expectation value of an observable M in an energy and J z - eigenstate \E : m) is defined as (M)j = 
(E, m\ M \ E,m) whereas for a TOP^-eigenstate ([ITT)) it reads 

WfioP. = TOP >(E,±\M\E,±f OP > 

= \ [{M) Jz + {TOP z MTOP z ) Jz ± (TOP z M>* 2 ± (MTOP z ) Jz \ . (13) 
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If M satisfies the anti-commutation relation {T0P Z1 M} = the first two terms of ([T3|) cancel each other and we 
obtain 

(M)± OP2 = ±\ [(TOP z M)\ - (TOP z M) Jm ] = T»Im ({TOP z M)J . (14) 

Since the expectation value of an observable is always a real number it follows 

{M)± OPz = 0, (15) 

i.e. the expectation value of such an operator M vanishes within any TOP^-eigenstate. 

In a similar manner one can calculate the expectation value of an observable TV, which commutes with J z and 
TOP z . Here expression (JT3|) yields 

(N)toP; = (N) Jz ±(NTOP z ) Jz . (16) 
Exploiting [iV, J z ] = and using {TOP Zl J z } = one obtains 

(NTOP z ) Jz = ^- (E, m\ NTOP z J z \E, m) = — — (E, m\ N J z TOP z \E, m) 

= -^(E,m\J z NTOP z \E,m) = -(NTOP z ) Jz (17) 
Hence (N TOP z ) Jz = leading to the final result 

(N)toP; = (N) Jz , (18) 

i.e. we obtain equal expectation values of N for a J z - or TOP^-eigenstate. 

Finally we provide some relevant expressions for the spatial probability density W in the different representations. 
The probability density belonging to the J^-eigenstate \E : m) = \u) ||) + \d) \[) reads 

W Jz (r,0,0) =r 2 sin<9 \\u(r,0)\ 2 + \d(r,0)\ 2 ] . (19) 

Since the whole ^-dependence is contained in phase-factors (see eq. (0) Wj z is invariant under rotations around the 
z-axis. For the TOP^-eigenstates ([ITT)) one obtains 

W T op 9 (rA<l>) = ! ^f^ [\u(r,6)\ 2 + \d(r,0)\ 2 + \u(r,n-6)\ 2 + \d(r,n-6)f (20) 
±2 {u (r, 0) d (r, tt - 0) + d (r, 0) u (r, tt - 0)} cos (2m 0)] . 

Wtop z (r, 0, (j)) is constructed from TOP z -eigenstates and obeys therefore the T0P 2 -symmetry. Since Wtop z (r, 0, <p) 
is a real scalar function the symmetries T and O are trivially satisfied leading to 

W TO p z = TOP z W TOPz (TOP*)- 1 = TP Z W TO p z (TP,)- 1 = P z W TO p z (P z )~\ (21) 
showing that Wtop z (r, 0, (j)) is invariant under z-parity. 



V. COMPUTATIONAL APPROACH 



We employ the linear variational principle [13 to find the eigenvalues and eigenvectors of the Hamiltonian (JJJJ. 
The idea is to expand the exact eigenfunctions of the Schrodinger equation in a complete basis set of spinor orbitals 
that converge efficiently and accurately towards the exact solution. Calculating the expansion coefficients results in 
a large-scale algebraic eigenvalue equation. 

The basis set I B^, Bz?\ we apply for given m takes the following form for the upper and lower component, 
respectively: 

(r,MI^) = i#'%)^r~W)ll> , {r,e,<p\B^=Rf k) (r)Y™ + H0A)\l) (22) 
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Here Y l m (0, <j)) denote the spherical harmonics. Due to the appearance of the upper and lower spinor orbitals for fixed 
m the energy eigenfunctions of the Schrodinger equation computed with this basis set are a priori eigenfunctions of 

j z . 

The radial part of the orbitals used for the expansion reads 



' fe)(r) = V (^fcj! e " ¥ (CO fc ^*(Cr) (23) 

with L^ k (r) being the Laguerre polynomials. The parameters k and ( can be tuned to optimize the convergence 
behavior in different parts of the spectrum. £ possesses the dimension of an inverse length. In order to converge 
eigenstates with the smallest possible number of basis functions l/( has to be adapted for every energy regime such, 

that it corresponds to the typical spatial extensions of the wavefunctions. For fixed k and ( the functions Rn' k \r) 
form a complete functional set in r-space but are non-orthogonal leading to a non-trivial overlap-matrix. The basis 
set (|22j> is complete in r-, 6- and spin-space. 

An eigenstate \E,m) of the stationary Schrodinger equation can now be expanded in terms of the basis-functions 

(G2) 

n<N,l<L h<N,l<L 

\E,m) = J2 a n,i\Bl)+ J2 b *j\ B ti) ( 24 ) 

n=0,Z=|m-i| n=0,[=|m+i | 

leading to the generalized spinor eigenvalue problem He = ESc, where H and S are the corresponding matrix 
representation of the Hamiltonian ® and the overlap matrix, respectively: 

oT I TJ\ oT \ / oT I TJ\ TdI 

H = I ) " f ( ) 7 ' f ( I S = I \ . x I (25) 





B l -\H\Bl lf ) (B l -\H\B 1 -) I I (B l ~j\B„ ll 

v n l i i n l J \ rd x 1 n'l' / J \ \ nl 1 n'Z 

The vector c contains the expansion coefficients of ((24)): 



6 fM 



(26) 



The two off-diagonal blocks of S vanish due to the orthogonality of the spin states | j) and ||). This does not hold for 
H since the magnetic field leads to a coupling of the two spinor components. With the basis set (|22|) all entries of 
the matrices ((25)) can be computed analytically. To do this we have exploited recurrence identities for the spherical 
harmonics and the Laguerre polynomials [13. The matrices H and S possess a particular sparse appearance, e.g. S 
is penta-diagonal, enabling us to go to large basis set dimensions. 

In order to solve the above generalized eigenvalue problem we employ the Arnoldi-method [l^ which is a so-called 
Krylov-space method. This approach is especially suited for eigenvalue equations involving large sparse matrices. The 
Arnoldi-Method is perfectly suited to provide low- lying excited states with fairly small basis sets. To calculate highly 
excited states without computing the lower excited states we use the so called shift- and-invert method. Here the 
generalized eigenvalue problem is transformed to (H — a S) _1 Sc f = \d . The eigenvalues E of the original eigenvalue 
equation are connected to A according to Ei = a + j. Thus eigenvalues lying close to the shift a become the 

largest eigenvalues in magnitude of the operator (H — aS) _1 S. Using the Arnoldi-method these eigenvalues lying 
in a predefined energy range of the spectrum will converge first. The shift- and- invert method allows in principle to 
compute arbitrary parts of the spectrum of the Hamiltonian (0) by resetting the shift a in successive calculations. 



Thereby the basis parameter £ has to be adapted. A good choice is £ = Q v 775 where E is the lower boundary of the 
corresponding energy range. We remark that there is no need to change the parameter /c, which we always put to 
zero. 

Due to the Hylleraas-Undheim theorem pjl] the numerical approximate eigenvalues provide an upper bound for the 
exact eigenvalues. In order to check the convergence behavior we calculate a number of eigenvalues Ef 1 for a given 
basis set size G\. In a second step the size of the basis set is increased significantly to G2 and the same eigenvalues 
are calculated again being now denoted by Ef 2 . As a measure of convergence we define the quantity 

rpGi fpG2 

Ki = &-Ehi (2T) 
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where the difference of the same eigenvalue for the two basis sizes G\ and G2 is divided by the distance to the lower 
neighboring eigenvalue. For Ki < 0.01 we assume the eigenvalue Ef 2 to be well converged. 

The maximum dimension of the basis sets we employed was approximately 17000. We thereby were able to converge 
several thousand eigenstates and eigenvalues up to energies corresponding to a hydrogen principle quantum number 
of n « 60 for weak gradients and gradients b ranging from 10 -10 to 10 covering the J^-quantum numbers me [— |, |] . 

VI. RESULTS 

In the following we discuss and analyze our results on the atom in the quadrupole field. We address for different 
regimes of the field gradient the effects of the quadrupole field on the spatial appearance of the wavefunctions as well 
as their spin properties. Furthermore we provide the selection rules for electric dipole transitions. A comparison with 
properties of the atom in a homogeneous magnetic field is performed to illuminate the specific features in a quadrupole 
field. 



A. Spectral Properties 

Depending on the gradient b of the quadrupole field one can distinguish three regimes: the weak, the intermediate 
and the strong gradient regime. In each of these regimes the atoms exhibits unique features. The notion of a 
weak/strong gradient does not refer to an absolute field gradient but depends on the degree of excitation of the atom. 
It is defined to be the regime, for which the magnetic compared to the Coulomb interaction is weak/strong. For 
weak gradients the level structure is dominated by the spin as well as the spatial Zeeman-term both of which depend 
linearly on the field gradient. The level degeneracies at b = are lifted and a linear splitting of the n-multiplets is 
observed. There is no overlap between neighboring n-multiplets which renders n an approximately good quantum 
number. Figure shows the linear splitting of the n = 20 multiplet for the states with quantum number m = 




FIG. 2: a: Energy level splitting of the m = ^-states of the n — 20-multiplet in the quadrupole field, b: Splitting of the same 
energies in a homogeneous field. 

The levels are arranged symmetrically around the energy for zero field. This is in contrast to the level splitting of an 
atom in a homogeneous magnetic field which is presented in figure|2t> for the same n- multiplet. Here the Zeeman-terms 
cause a linear energetical splitting into only two branches. The latter (for fixed m) is the result of the two possible 
spin orientations in a homogeneous magnetic field. In a quadrupole field L z and S z are not conserved separately and 
an assignment of spin orientations is therefore not possible. 

With increasing gradient the diamagnetic term with its quadratic dependence on b (or B in the homogeneous field) 
becomes increasingly important thereby resulting in a non-linear behavior of the energy curves. The intra n-manifold 
mixing regime where different n-multiplets are still energetically well separated but different angular momentum states 
(/-states) mix can be shown to scale as b cx n -6 . In the high gradient regime different n-multiplets begin to overlap 
and avoided crossings dominate the spectra. The onset of this inter n-manifold mixing scales according to b cx n~~ 
whereas in the homogeneous field the scaling is B cx n~2. Figure 03 shows the energy levels of n-multiplets ranging 
from n = 33 to n = 37 with m = ±| in the regime < b < 10 -8 . For such a degree of excitation the intra n-manifold 
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FIG. 3: Energy spectrum of the multiplets n — 33 — 37 for the states with m = ±|. For low gradients the energy levels split 
almost linearly followed by a transition region where the diamagnetic term becomes increasingly more important. For high 
gradients different n-multiplets overlap and no symmetries, i.e. approximate quantum numbers, are left. 

mixing sets in at b « 5 • 10 -10 . Therefore the linear splitting due to the Zeeman-term as shown in figure is hardly 
visible. 



B. Low lying states and large gradients 

Although gradients b ^> 10 -10 are not experimentally accessible it is at this point instructive to illuminate their 
impact on low-lying electronic states. This way some principle features of the influence of the quadrupole field on the 
electronic structure of the atom will become clear. 

Here we consider the behavior of the lowest three n-multiplets as a function of b for < b < 10 -2 (see figure 0J. 
The splitting of the n = 3 multiplet is linear up to b = 5 • 10 -4 , where inter n-manifold mixing begins. 

The kink of the energetically highest-lying curve atfr = 2-10 -3 is the result of an avoided crossing with an energy 
level from the upper n = 4 multiplet. Moving to higher gradients we encounter further avoided crossing. Crossings 
between levels belonging to different quantum numbers m are exact, since the corresponding states belong to separate 
symmetry sub-spaces. The n = 2 multiplet does not show any indication of inter n-manifold mixing throughout the 
whole gradient range considered here. Figure 0| shows also that the energy of the degenerate ground state depends 
very weakly on the gradient b: We observe a maximum relative energy shift of about 0.04% at b = 10 -2 . 
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FIG. 4: Energy level structure of the three lowest n-multiplets. The J z quantum number ranges from 



to 



C. Radial compression of the wavefunctions 



The operator r obeys the commutators [r, TOP z ] = and [r, J z ] = 0. Together with (p"8j) this yields 

(r) = (r)rop z = ( r )j z • 



(28) 



Thus the expectation value of r in the J z eigenstates is the same as in the T0P 2 -eigenstates. 

For the one-electron problem without external field we have (r) H = \ (3n 2 — /(/ + 1)). Since I satisfies < I < n — 1 

one obtains corresponding upper and lower bounds for the expectation value of r: n (n + |) < (r) H < §n 2 . In figure 
[5] this expectation value is shown for the atom in the quadrupole field as a function of the principle quantum number 
n for two different gradients (n serves here as a label for the energy levels and is not a good quantum number!). 
The radial expectation values lie in between the boundaries given by the field-free inequality indicated by the dashed 
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FIG. 5: Expectation value of the radial coordinate r plotted against the principle quantum number n at different gradients (a: 
b= 1CT 10 , b: b= 10" 8 ). 

lines (see figure Et). Expectation values of states belonging to the same n-multiplet are arranged in vertical lines 
expressing their energetical degeneracy (see figure Et for b = 10 -10 ). Here states with large expectation values of r 
possess small expectation values of the angular momentum and vice versa. The situation is different for b = 10 -8 . 
Here a systematic decrease of the r-expectation values takes place for states with an energy corresponding to n > 30. 
This is also the energy regime where the inter n-manifold mixing sets in. Here the mentioned pattern of vertical 
lines dissolves and a irregular looking distribution emerges. Although the expectation values are not located between 
the dashed boundaries the distribution is still widespread which can be ascribed to the different angular momentum 
distribution of the states. 

We now investigate the deformation of the ground state wavefunction in both the quadrupole and the homogeneous 
field. Due to their simple nodal structure they are suited best to point out differences arising for the two different field 
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configurations. To demonstrate the deformation effects on the electronic cloud we consider for reasons of illustration 
an extremely high gradient. Figure shows the spatial probability density of the ground state in the quadrupole field 




r [a. u.] r [a. u.] 

FIG. 6: Spatial probability density of the ground state in the quadrupole field (a) and a homogeneous field pointing in the 
z-direction (b) for b = B = 10 and m = \. 

JSt) at b = 10 and the homogeneous field ^jp) at B = 10 (m = - in both cases). 

For the quadrupole field we observe an asymmetric deformation with respect to the = ^ -plane: the electronic 
wavefunction is almost completely confined to the upper half- volume (0 < §) which is a consequence of the symmetries 
discussed in section ITTT1 z-parity is not conserved and eigenstates appear in pairs one being the mirror image of the 
other with respect to the reflection at the x — y-plane. Furthermore we observe that the electronic motion is localized 
particularly along the two 'channels' for = corresponding to the lower z-axis and 6 = ^ being the x — y-plane. 
This property as well as the detailed shape of the electronic probability density in the individual half-volume is 
determined by the diamagnetic term which is dominant in the high gradient regime. For the quadrupole field it is 
proportional to sin 2 cos 2 reaching its maximum value at = ~ , ^ . The probability density in a homogeneous 
field (see figure EJ)) exhibits the above-mentioned corresponding reflection symmetry due to the invariance of the 
corresponding Hamiltonian with respect to z-parity. Here we observe the maximum of the probability density at 
= ^ and a deformation towards = tt and = leading to a cigar-like shape. The diamagnetic term is proportional 
to sin 2 having its maximum value at = ^ thus coinciding with the regions possessing the strongest deformation 
of the probability density. For both field configurations the probability density vanishes at r = 0. The maximum 
is reached at r « 0.75 in the quadrupole field and for r ~ 0.5 in the homogeneous field followed by an exponential 
drop-off. 



D. Properties of the electronic spin 

1. Expectation values 

In a homogeneous magnetic field the projection of the spin operator onto the direction of the magnetic field is a 
conserved quantity. In this case one can choose the energy eigenstates to be also eigenfunctions with respect to S z . 
This restricts (S z ) to the two values ±~. In the quadrupole field S z is not conserved, and we therefore consider the 
expectation value (S z ) of the electronic states. We have 

(S z ) Jz = \[(u\u)-(d\d)}. (29) 

Figure Q shows the distribution of {S z )j z for electronic states of the m = \ subspace as a function of the principle 
quantum number n. Since S z is not conserved the values of (S z ) are allowed to cover the complete interval [— ^, |] . 
For b = 10 -10 (figure E^) an d a low degree of excitation the expectation values are evenly distributed over the interval. 
When reaching highly excited states this pattern becomes increasingly distorted. The expectation values agglomerate 
at —0.35 and 0.35 for n > 50. Due to the approximate degeneracy of the energy levels at low energies, i.e. small n, the 
values of (S z ) form vertical lines. These lines widen for higher n. Since at b = 10 -10 no significant n-mixing up to our 
maximum converged energy levels takes place, neighboring lines are well separated. For a higher gradient b = 10 -8 
(figuretZb) the above properties are equally present for low-lying states. However, with increasing excitation energy we 
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FIG. 7: Expectation values of the z-component of the spin operator as a function of the quantum number n for different 
gradients (a: b = 10" 10 , b: b = 10" 8 ). 

now observe a complete n-mixing regime where the regular patterns disappear and we obtain an irregular distribution 
of (S z ). Overall the distribution narrows, e.g. for n = 40 the occupied interval is approximately [—0.3,0.3]. 

We remark that due to the fact that the ^-operator anti-commutes with TOP z we have (S z )^ OP = 0. Apparently 
there is no preferred direction for the electronic spin in a state obeying the TOP z symmetry. 

2. Spin polarization 

Due to the coupling of the spatial to the spin degrees of freedom we have to expect a dependence of the spin 
orientation on the spatial coordinates. To investigate this in more detail we study the spatially dependent S z - 
polarization Ws(r). For a J^-eigenstate \E,m) it reads 



W s {f) = 



(E,m\r)S z (r\E,m) _ 1 \(u | r)\ 2 - \(d | r)\ 2 



(E,m\r)(r\E,m) 2 \( u | f)| 2 + \( d \ f)| 



2 • 



(30) 



Figure [SI shows the spatial probability distribution (J^t) and the ^-polarization JSb) for the 45th excited state for 




50 r ,100 
r[a.u.] 



50 r ,100 
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FIG. 8: a: Spatial probability density of the 45th excited state for m = \ and b — 10 8 in logarithmic representation, b: 
S z -polarization for the same state. 



m = | and b = 10 -8 emerging from the n = 7 multiplet at b = 0. The probability density is shown for a logarithmic 
scale to provide details for small values. In contrast to the constant ^-polarization we would encounter in the absence 
of a field or a homogeneous field we observe a complex pattern of domains exhibiting different spin orientation (white: 
spin up, black: spin down). At low r values these domains form a pattern similar to that of a chess board. With 
increasing r a transition region is encountered where the formation of stripes with different spin orientation begins. 
The junctions where four spin domains meet each other coincide with the nodes of the spatial probability density. The 
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Coulomb interaction as well as the spin-Zeeman term are responsible for the interwoven network of island of different 
spin orientation. The additional presence of the orbital Zeeman and the diamagnetic term leads to a deformation of 
this network. 

Highly excited states, i.e. Rydberg states, typically possess a large spatial extension. In figure EH the spatial 
probability density an d the ^-polarization (03) are depicted for the 1117th excited state of the m = | subspace 
(b = 10 -8 ). At small r- values the 5^ -polarization shows the chess board pattern, we have already observed for the 
low-lying excited states, followed by a striped region. Whereas for low radii the spin orientation changes locally from 




FIG. 9: Spatial probability density (a) and -polarization (b) for the 1117th excited state for m — \ and b = 10 8 . At large 
r the SVpolarization becomes similar to Wg indicating an antiparallel alignment of the electronic spin to the magnetic field. 



island to island generating an appealing pattern we observe an overall tendency of the electronic spin polarization 
in the region characterized by stripes (r « 2000). Here, independently of the nodal structure, the spin orientation 
changes smoothly from downwards at = tt to upwards at = 0. This feature can be understood by inspecting the 
spin Zeeman term only. The corresponding Hamiltonian reads 

H S = -^(r) = -r^ s[n0ei , 2cQs0 y (31) 

Here the complete dynamics takes place in spin space since the spatial coordinates (r, 0, <p) are entering as parameters. 
The eigenvalue problem belonging to ((31)) possesses the two solutions 

/ l cos cos^ t) „—id> \ 

$±(r,0,(t>)= i ™l J (32) 

with the energies 

E± = =F^r\/l + 3cos 2 l9 = ^\fl\\B(r)\. (33) 



These energies correspond to those of a spin oriented parallel (E-) or antiparallel (£"+) to the magnetic quadrupole 

7 S 



field. Constructing the ^-polarization W^r of the eigenstates (|33[) yields 



_ 1 (2 cos (9 ± Vl + 3 cos 2 (9) 2 - sin 2 
5 ~ 2 (2 cos ± Vl + 3 cos 2 (9) 2 +sin 2 (9' 



W£ = ^ . . u 5 o.o — T7- (34) 



Both W£ and Wg are shown in figure They do not depend on the radial coordinate and the azimuthal angle 
The ^-component of the spin is oriented upwards/downwards at 9 = and points in the opposite direction at = 



For 6 = 7^ there is a transition region with W s being close to zero. 

This is precisely the above-discussed behavior we observe at large r for the ^-polarization of the Rydberg state 
depicted in figure |2t>. Apparently in this special state the electron spin prefers an antiparallel alignment to the 
external field at large radii which corresponds to Wg. We remark that similar results on the spatial spin polarization 
are obtained if one considers the quantity SB(f) instead of S z . 
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FIG. 10: S z -polarization of the two eigenstates of the Hamiltonian (|31|) . The electronic spin either points parallel (W s ) or 
antiparallel (Wg) in the local direction of the field. 

E. Electromagnetic transitions 

The interaction of an atom in the magnetic field with external electromagnetic radiation leads to transitions among 
its electronic states which we describe employing the dipole approximation. The amplitude for a transition from 

the initial state \i) to the final state |/) is given by the matrix element ,7V = (i\D a ±^ \f). Considering ir and 
<j ± transitions the operator takes the common forms D n = z = rcosO and D a ± = (x ± iy) = -^rsmOe^ 1 ^, 
respectively. The symmetries of the J z and TOP z eigenstates result in selection rules which we will discuss in the 
following. 

Evaluating the transition matrix element for 7r-transitions (E f , m'\ r cos 6\E,m) leads to the selection rule m! — m = 
0. For <j ± -transitions the corresponding matrix element reads (£", m'\ r sin# e ±% ^ \E, m) which is only non-zero if 
m! — m — ±1. Hence, electromagnetic transition between J z states occur only if Am = 0, ±1 reminiscent of the 
situation in a homogeneous or even no magnetic field. The corresponding matrix element for 7r-transitions between 
TOP^-eigenstates is (£", H'\ r cos 9 \E, Ft) leading to the selection rule Ft 7 ^ Ft, i.e. only transitions between states 
with opposite TOP z symmetry are allowed. There are no selection rules for <j ± -transitions between TOP^-eigenstates. 
The selection rules are summarized in table [I] 



transition type 


J z — states 


TOP z - states 


Am — m' — m 


An = IT - n 


7T 





±2 


a+ 


1 




<J~ 


-1 





TABFE F Selection rules for dipole transitions between J z -eigenstates and between TOP^-eigenstates. 

We have calculated the dipole strengths for transitions from the ground state (m = |) to excited states with m = | 
(^-transitions (figure ITTjl ) and m = | (a~ -transitions ( figure IT2jl). For comparison we also present the transitions 
for an atom in a homogeneous magnetic field indicated by dashed lines whereas transitions in the quadrupole field 
are represented as solid lines. For a detailed discussion of electromagnetic transitions in the homogeneous field we 
refer the reader to @, 113 EH- I n order to make both results comparable we have chosen the gradient and the 
magnetic field strength such that we observe approximately equal splitting of the n-multiplets with increasing energy. 
Specifically we chose b = 10~ 7 and B = 10 -4 . Figure ITTk shows the dipole strengths for 7r-transitions. Here the 
line possessing the largest wavelength emerges from the transition to the n = 10 multiplet which already exhibits 
significant features of intra n- manifold mixing. From here up to a certain wavelength (« 91.5nm) neighboring lines 
are well separated. Transitions with smaller wavelengths involve levels showing inter n-manifold mixing becoming 
visible by the overlapping of neighboring groups of lines. In figure ITTb we provide a higher resolution picture of 
the transitions to the n = 12-manifold in the intra n-mixing regime. Apparently in both cases each main line is 
accompanied by a series of sub-lines. In the quadrupole field besides the main line at A « 91.803nra there are four 
major sub-lines situated at the outer edge A ~ 91.81nm. We observe a number of sub-lines almost equal in height 
compared to the main line for the homogeneous field. We note that the sub-line possessing the maximum strength 
always belongs to a transition in the quadrupole field. An overall feature is the fact that a significantly lower number 
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FIG. 11: Dipole strengths of Am = tt— transitions from the ground state to excited levels. Solid lines denote transitions in 
the quadrupole field whereas transitions in the homogeneous field are indicated by dashed lines, a: The line possessing the 
largest wavelength belongs to the n — 1 — > 10 transition, b: Magnification of the line structure of the n = 1 — > 12 transitions. 
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FIG. 12: Dipole strengths of Am = — 1 a— transitions (a "-transitions) from the ground state to excited levels. Solid lines 
denote transitions in the quadrupole field whereas transitions in the homogeneous field are indicated by dashed lines, a: The 
line possessing the lowest wavelength belongs to the n — 1 — > 10 transition, b: Magnification of the line structure of the 
n — 1 — > 12 transition. 



of sub-lines occurs in the homogeneous field compared to the quadrupole field: Symmetry properties deriving from the 
conserved z-parity give rise to additional selection rules and lead to a reduction of the number of allowed transitions. 
In figure [T21 we present the transitions to states m = | starting from the transition to the n = 10 multiplet. Again 
the change from the inter to the intra n- mixing regime is observed (figure Et)- Here the threshold wavelength at 
which the overlapping of adjacent groups of lines starts is also at A « 91.5nra. Compared to the wavelengths in 
the homogeneous field we find the transitions in the quadrupole field to be systematically shifted towards larger 
wavelengths. In figure fT21) we show the magnification of the line structure for the n = 1 —> 12 transition. In the 
homogeneous field each transition to a fixed n-multiplet is dominated by a major line accompanied by a group of 
significantly weaker sub-lines. In contrast to this we find many transitions of equal strength in the quadrupole field. 
This is in contrast to the above-discussed 7r-transitions. The total number of transitions is again much larger in the 
quadrupole field compared to the homogeneous field. 



F. Magnetic field induced permanent electric dipoles 

In a homogeneous magnetic field parity is a symmetry and therefore atomic electronic eigenstates (for fixed nucleus) 
do not possess a permanent electric dipole moment. Let us now investigate the electric dipole moment of the electronic 
states in the quadrupole field. Due to the selection rule given above the expectation value of D a ± in the J^-eigenstates 
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is 

(D a ± ) = (E, m\ r sin 6 e ±icf) \E,m) = 0. (35) 
However, the expectation value of is in general non-zero. (D^) is shown in figure IT^T a. b) for the two gradients 




n n 

FIG. 13: Expectation value of the dipole operator plotted against the principle quantum number n for different gradients 
(a: b= 1(T 10 , b: b= 1(T 8 ). 

b = 10 -10 and b = 10 -8 , respectively. For b = 10 -10 the electric dipole moments belonging to the same n-multiplet 
are arranged along vertical lines, which is a result of the approximate degeneracy of the energy levels. The substates 
for fixed m of a given n-multiplet exhibit different dipole momenta spreading between two (upper and lower) bounds 
that depend linearly depending on n. States with a large electric dipole moment emerge from field-free states (with 
increasing b) that possess small values for the angular momentum and vice versa. With increasing gradient and degree 
of excitation the n-mixing starts and disturbs the observed regular pattern. For b = 10 -8 and n > 35 the distribution 
of the dipole moments becomes completely irregular. Therefore we encounter the remarkable effect that the external 
magnetic quadrupole field induces a state dependent permanent electric dipole moment. This is the result of the 
asymmetric form of the wavefunction in the quadrupole field (see discussion in section IVl CI and figure EJ. Exciting 
the atom from the ground to a corresponding excited state via lasers with the corresponding transition frequency it 
is therefore possible to prepare an atom with a desired permanent electric dipole moment ! 

VII. CONCLUSION AND OUTLOOK 

We have presented a detailed study of the electronic structure of atoms in a magnetic quadrupole field. Opposite 
to the common adiabatic description of the atomic center of mass motion in inhomogeneous magnetic fields we have 
focussed on the internal electronic states of the atoms exposed to the inhomogeneous field. We have employed an 
effective one-particle approach including both the coupling of the electric charge and the magnetic moment (spin) to 
the field. A spinor-orbital based method to compute the eigenfunctions of the stationary Schrodinger equation has 
been developed and applied. We have utilized a 'Sturmian' basis set to study several thousands of excited states 
for a regime of gradients of more than 10 orders of magnitude. Due to the inhomogeneity of the quadrupole field 
the spatial and spin degrees of freedom are coupled in a unique way. As a result the system is invariant under a 
number of symmetry operations acting on both degrees of freedom. These are the unitary symmetries related to 
the conservation of the total angular momentum J z and the discrete operation P$OP z as well as the anti-unitary 
generalized time reversal symmetries TOP z and TP^. These operations constitute a non-Abelian symmetry group 
which leads to a two-fold degeneracy of each energy level. This is a remarkable feature in particular since it occurs in 
the presence of the external field. Without a field such (Kramers) degeneracies of a spin- ^-particle are well-known. 

We have shown energy spectra up to excitation energies corresponding to principal quantum numbers of n « 60. 
The spectrum exhibits distinct characteristics for the weak, intermediate as well as the strong gradient regime. At 
weak gradients a linear splitting of the energy levels is observed. For given m the sublevels of an n-multiplet split 
symmetrically around the zero field energy. For the homogeneous field a splitting in only two branches is observed due 
to the two energetically different orientations of the electronic spin. Employing higher gradients adjacent n-manifold 
are still separated but we encounter intra n-manifold mixing. Scaling relations for the onset of the intra- as well as 
the inter n-manifold mixing in a quadrupole field have been given. To understand some general features of the effects 
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generated by the quadrupole field we have studied the spectrum of low excitations in very strong (experimentally 
not accessible) gradients. We have found that the ground state energy virtually stays the same up to gradients of 
b = 10 -2 whereas the first few excited states already experience severe changes. The latter obviously only holds 
for hydrogen. However, one way to proceed is to investigate Rydberg states of nonhydrogenic atoms in detail. Here 
especially scattering with the inner electron shell leads to modifications of the energy spectra which can be understood 
by quantum defect theory 0, [2^ . At least for Rydberg states which possess a large angular momentum we do not 
expect significant modifications induced by core scattering processes. 

In our investigation of the electronic spin properties we found the ^-expectation values of the electronic states to 
form a regular pattern at low gradients. For higher gradients and/or higher excitations we observe a transition to an 
irregular, much narrower, distribution. In order to analyze the spatial behavior of the electronic spin we introduced 
the 5^ -polarization. Due to the unique coupling of the spatial and spin degrees of freedom this quantity reveals a 
rich nodal structure. At low radii we observe a chess-board like pattern of islands with alternating spin polarization 
changing to a striped pattern at large radii. For Rydberg states the chess board pattern changes to an irregular 
network of islands exhibiting opposite spin polarizations. Within the striped region an envelope behaviour of the spin 
orientation on the angle is observed. To investigate this phenomenon we have studied a Hamiltonian describing 
the coupling of the electronic spin to the magnetic field, neglecting the spatial dynamics. Analytically calculating its 
eigenstates we could reproduce the above mentioned ^-dependence of the 5^ -polarization of Rydberg- states for large 
radii. 

We have derived selection rules for electromagnetic dipole transitions among eigenstates to the J z as well as TOP z 
operator. Linearly polarized transitions take place only between states with the same quantum number m or states 
with different TOP z symmetry. Furthermore we have calculated the amplitudes for tt- and o~ -transitions emerging 
from the ground states (m = |) to levels lying in the intra and inter n-manifold mixing regime. We have compared 
the results with those of a homogeneous magnetic field discovering significant differences in the line structure. Due 
to additional selection rules arising from the conservation of z-parity fewer transitions are observed in a homogeneous 
magnetic field. Calculating the expectation value of the dipole operator it turned out that the electric dipole moment 
of the electronic states in the quadrupole field is in general nonzero. For low degrees of excitation and low gradients 
we observed an almost linear increase of the maximum dipole moment of the n-multiplets. This pattern becomes 
increasingly distorted when moving to higher gradients and/or a higher degree of excitation. The no n- vanishing 
dipole moment is a result of the spatially asymmetric arrangement of the electronic wavefunction with respect to the 
x — y-plane. 

In order to achieve the experimental realization of the discussed system the so called atomchip 6] seems to be 
a promising device. Here current carrying nano- structures with dimensions in the /im-regime can generate high 
gradient quadrupole fields. At the moment gradients such as b = 10 -10 or even slightly above represent certainly 
the limit. However, almost all effects discussed here are not due to the diamagnetic interactions but have their 
origin in the interplay between the Coulomb as well as the spin- spatial Zeeman terms. Therefore weak gradients 
should not represent a principle obstacle to experimentally observe the derived properties. In particular the magnetic 
field-induced permanent electric dipole moments would potentially find applications in e.g. quantum information 
processing |2g, |2jJ |2S ui| • Exciting the atom from the ground to a corresponding excited state via lasers with the 
corresponding transition frequency it is possible to prepare an atom with a desired permanent electric dipole moment. 
In order to build a working two qubit gate one is interested in finding ways to establish a controlled interaction 
between two qubits to gain a phase shift. Considering a trapped Rydberg state as a qubit the interaction between 
two of them can be realized via dipole-dipole interaction. This interaction could be switched on and off on demand 
by changing the atomic state in the quadrupole field. 

Taking into account the finite mass of the nucleus it has been shown ^J^J^^ that the center of mass (CM) and 
electronic motion do not separate i.e. do not decouple in the presence of a homogeneous magnetic field. To enter 
the corresponding regime where the residual coupling becomes important certain parameter values (excitation energy, 
CM energy etc.) have to be addressed. A variety of intriguing phenomena due to the mixing of the electronic and 
CM motion have been observed consequently. Examples are the classical diffusion of the CM due to its coupling to 
the chaotic internal moti on 1 131. the giant dipole states of moving atoms in magnetic fields [l^ as well as the the 
self-ionization process l23ll24j due to energy flow between the CM and electronic degrees of freedom. In the present 
investigation we have fixed the position of the CM of the atom at the center of the quadrupole field i.e. we have 
adopted an infinitely heavy nucleus. On the other hand assuming that an experimentally prepared atom is ultracold 
will certainly minimize the CM motional effects. Nevertheless, a residual coupling is unavoidable and its impact on 
the electronic structure is, at this point, simply unknown: A full treatment of the two-body system certainly goes 
beyond the scope of the present investigation and requires both from the conceptual as well as computational point of 
view major investigations. However, one should note that the symmetries discussed here equally hold for the moving 
two-body system i.e. the total angular momentum is conserved and the unitary as well as antiunitary spin-spatial 
symmetries, now applied to both particles, are also present. 
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